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Abstract. We investigate the nature of the plasma analogy for the Laughlin wave function on a torus 
describing the quantum Hall plateau at We first establish, as expected, that the plasma is screening 

if there are no short nontrivial paths around the torus. We also find that when one of the handles has a short 
circumference — i.e. the thin-torus limit - the plasma no longer screens. To quantify this we compute the 
normalization of the Laughlin state, both numerically and analytically. In the thin torus limit, the analytical 
form of the normalization simplify and we can reconstruct the normalization and analytically extend it back 
into the 2D regime. We find that there are geometry dependent corrections to the normalization, and this 
in turn implies that the plasma in the plasma analogy is not screening when in the thin torus limit. Despite 
the breaking of the plasma analogy in this limit, the analytical approximation is still a good description of 
the normalization for all tori, and also allows us to compute hall viscosity at intermediate thickness. 
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1. Introduction 

The Laughlin wave function is the drosophila of representative trial wave functions for the fractional quantum 
Hall effect (FQHE). It was introduced in a planar geometry by Laughlin [Lau83| . generalized to the sphere 
by Haldane |Hal83a| . to the torus by Haldane and RezayijHRSSj, and has has been extensively studied 
on these geometries ever since. Most remarkably, its elementary excitations - its quasi-particles - have 
fractional charges |Lau83| . and are believed to obey anyonic statistics |ASW84] . The latter, although not 
proven experimentally, is strongly suggested by extensive numerical studies |KM99I lKL99| . Recently, matrix 
product states have been successfully used to probe the anyonic statistics [ZMP 13l 1C V 13| . 

The theoretical argument for fractional statistics is based on the so called plasma analogy, first introduced 
by Laughlin |Lau83| . The main observation is that, introducing a particular factor dependent on the quasi¬ 
particle positions into the normalization of the wave functions, the corresponding full normalization constant 
can be expressed as the partition function for a classical screening plasma. In this analogy the quasi-particles, 
at positions rji, appear as test charges with a strength one of the electron charge. As the plasma is 
screening, the partition function does not depend on the positions of the quasi-particles, provided they are 
separated much further than the screening length. From this, one concludes that the full ? 7 i-dependence of 
the normalization is accounted for, and from this it is straightforward to calculate the fractional statistics 
parameter. 

Another way of saying this, is that the "holonomy" i.e. the phase related to adiabatically dragging one 
quasi-particle at rji around another at 772 will equal the "monodromy", i.e. the explicit phase obtained from 
the normalization factor. For this to occur, no additional contribution from any Berry phase related to the 
adiabatic change must be accumulated. 

Since screening is essential for the argument of vanishing extra Berry phases, we may ask the contrary 
question: When does the plasma not screen? Can the screening properties of the plasma change if the quasi¬ 
particles move too close to the edge of the sample or if there are non-trivial loops on the spatial surface with 
short length scales |KM991 IKL99] ? A related question that is also important for the Laughlin states is the 
following: How does the free energy of the plasma depend on the topology, geometry and size of the surface 
on which the QH liquid resides? This is particularly interesting in the case of a torus where the normalization 
can depend on the aspect ratio of the two principal axes, t. t is a complex valued parameter that can be 
changed while still keeping the density and number of electrons fixed. In the context of the plasma analogy, 
this corresponds to asking if the plasma is screening for all r. 

Recently the r-dependence of the normalized Laughlin wave function, ipL', has been determined using 
techniques from conformal field theory(Rea09]. Later numerical work |RR 11| confirmed that there where no 
extra r-dependent contributions to this normalization. The numerical analysis was however only performed 
for special values of r, so the question remains if there is need for corrections when considering general r. 

It is fruitful to think of t as a parameter that can be adiabatically changed, in the same way as the 
quasi-particle positions can be adiabatically dragged around each other. Under such an adiabatic change 
analogous questions appear as for the quasi-particles, i.e. is the Berry phase the same as the monodromy 
of changing r. In fact, for the Laughlin state this computation can be considered even more interesting as 
changing r can transform the q degenerate Laughlin states into each other, giving rise to Berry matrices. 
These matrices, which are well defined even without the presence of quasi-particles, would form a non-abelian 
representation of the modular group, the mapping class group of the torus, while exchanges of quasi-holes 
only give rise to phase factors. 

In this paper we study the normalization properties of with no quasi-particles present. For general 
r we obtain the normalization as a finite - but intractably large - sum in powers of We expand "0^ 

in T by rewriting it in a basis of single particle orbitals i.e. a Fock basis. Progress has been made on the 
plane and sphere by identifying the Laughlin states with Jack polynomials jBHOS]. The Jack polynomials 
can be recursively calculatedjLLMQQ] by starting with the root partition in the orbital occupation basis and 
applying squeezing rules. The root partition is the configuration of electrons where the distance between all 
occupied orbitals is maximized. For the Laughlin state at = |, the root partition is the orbital occupation 
pattern ... 1001001001.... 

On the torus, no known squeezing rules exist and the Laughlin state can not be written in terms of 
Jack polynomials. On the other hand, there is a well defined limit, the Tau-Thouless (TT) limitjBKOS], in 
which the Laughlin state is precisely a single slater determinant corresponding to the root partition. From 
the TT-limit it is possible to generate the Laughlin state at a general aspect ratio r using a differential 
equation in this parameter jZNSlS] . The differential operator acts on the Fock expansion of the Laughlin 
state. In the TT-limit the analytical expression for the normalization also simplifies greatly can be expanded 
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Figure 1. The relationship between the Cartesian coordinates (x^y) and the dimensionless coordinates 

(a:, y). In the figure one can also see that ri = is interpreted as the skewness, and T 2 = ^ as the aspect 

L'x 

ratio, of the torus. The area of the torus is fixed to be L^Ly = L^T 2 = 2iTN^tg. 


perturbatively in powers of 

To find the Fock expansion of the wave function is useful since it would allow for a direct computation 
of the normalization of the Laughlin state as well as many other interesting quantities! It is conjectured that 
for a torus large enough, the r-dependence of the Laughlin state is fully captured by the normalization factor 
proposed by Read |Rea09| . In that case, there would be no extra dependence on the precise geometry of the 
torus, just as the normalization is insensitive to the positions of well separated quasi-holes of the Laughlin 
state |Lau83| . if the factors mentioned before are taken into account. 

For small system sizes the normalization of the Laughlin state can be computed exactly for any value of 
r. The rough picture is that when the torus is large enough in both directions then there is no r-dependence 
in the normalization. However, if one of the torus handles becomes too thin there is r-dependence in the 
normalization. See e.g. figure ([^ in Section This is consistent with the intuitive picture that at some 
length scale the screening properties of the Laughlin plasma analogy break down. 

In this paper we first analytically extract the Fock coefficients for the Laughlin state in a generic 
torus geometry r. We obtain a sum that is finite, but untractably large for practical purposes. From 
the analytical expressions we approximate the leading behavior of the Fock coefficients in the TT-limit, 
and using this leading order expansion we comment on the validity of the plasma analogy in the TT-limit 
and on the universality of the Hall viscosity. We find that the plasma analogy is no longer screening when 
the circumference of the torus becomes sufficiently small, and we compute the leading corrections to the 
normalization in this limit. We conclude that although the plasma is screening in the thermodynamic limit, 
it will always cease to to screen if there is a loop around the torus shorter than some characteristic length 
scale. Intuitively, this distance, which we find to be roughly 6 magnetic lengths, should be on par with the 
screening length of the plasma. 

This paper is organized as follows. Section [^contains a brief summary of relevant notation. In Section]^ 
we introduce the analytical Fock expansion for the Laughlin state on the torus. In Section]^ we use the Fock 
expansion to study the Laughlin state in the TT-limit. We address both the proper norm in the TT-limit as 
well as extrapolate back to the thick torus. We also show that the quantum Hall viscosity of the Laughlin 
state is different in the TT-limit than in the thermodynamic limit. Finally, in Sectionl^we sketch the route to 
obtain a Fock expansion for the chiral Haldane-Halperin hierarchy states |Hal83al IHELl83b| constructed using 
CFT techniaues jFHSli] . 

2. Lowest Landau Level Wave Functions 

For self consistency, but also to define our notation, we give the basic formula for the lowest Landau level 
(LLL) wave functions and translation operators on the torus. 

In this paper we will exclusively work in the dimensionless coordinates where 0 < a;, y < 1, defined on a 
unit square. They are related to the dimensionful coordinate z as z = L {x + ry). See Figure]^ The complex 
parameter t = ti + it 2 parametrizes the torus geometry and is defined in the complex upper half plane. 

The single particle Hamiltonian is the Landau Hamiltonian 
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In these units the vector potential is {Ax,Ay) = (2TTN^i'^By,0^ = (T2L^By,0y We have here chosen to 
use, what we call the r-gauge, as in the physical coordinates z = x + ly the vector potential is perpendicular 
to T = (ri,T 2 ). The physical area of the torus is A = T 2 L^ = 2ttN^£‘^q and is penetrated by magnetic 
fluxes; we also deflne e = ^. We will throughout this paper let the magnetic length be = 1. The Landau 

Hamiltonian can be brought to the form of a harmonic oscillator H = huj (a^a + 5 ) by introducing ladder 
operators 


a = ^02 + ^Ly'j 

a'l' = - V2 (^^2 - ^Ly'^ ■ 


The energy levels are called Landau levels (LL) and are all iV^-fold degenerate. Acting with a on a state in 
the lowest Landau level (LLL) gives the equation (^2 + V’lll = 0. The solution to this equation shows 
that all LLL wave functions have the form 


f’LLL = f (z), 


where f (z) is a holomorphic function. 

The translation operators that commute with the Hamiltonian are finite magnetic translation operators. 
There are two minimal translation operators tx and ty, that move a state a finite distance along the principal 
axes of the torus, i.e. in the directions x and y. The distance translated by tx is e and so tx* makes a full 
revolution around the torus. The periodic boundary conditions are then defined through the equation 



where the (fa can be thought of as fluxes threading the two handles of the torus. The translation operators 
expressed in the r-gauge are tx = and ty = From these two minimal translations we construct 

the full set of translation operators as 

j. _ emdx-\-€.ndy-\-i2'Knx 

^ 


The minimal translations are identified as tx = ti,o, ty = to.i- The operator tm,n moves a wave function 

a distance z ^ z + eL fm + rn) and has commutation relation tm,ntm',n' = ^^tm’,n'tm,n- 

We choose to diagonalize the single particle wave functions in the LLL with respect to tx- For periodic 
boundary conditions {(fx = <fy = the single particle basis states are 


Vk (z) 




ek 

0 





( 1 ) 



nrTN^(y+t+ek)‘^ 


with properties txTjk = and tyPk = Pk+i- The eigenstates for generic boundary conditions are trivially 

obtained by acting with t_^y on ijk- There are Nif, single particle orbitals, one for each flux quantum. 
Since the LLL has a non-commutative geometry there is a connection between the momentum k and the 
expectation value of the position in y-direction. This relation is simply (y) = —ek. As such, the momentum 
label can alternatively be thought of as a physical displacement label of a one-dimensional system. This view 
is especially fruitful in the TT-limit where the width of the orbitals is much smaller than the inter-orbital 


distance, ay <C e. Here ay = \J (y^) — (y)^ is the standard deviation measure of uncertainty in the y-direction. 
The i?-function appearing in the above formula is defined as 




a 

b 


{z\t) 


00 

'y ^iTrT{t+a)^ ^z2Tr{t+a){z+b) 


t——oo 


( 2 ) 


and is frequently encountered when considering wave functions on the torus. 

Note that, in Q due to the holomorphic structure of ZllLj the exponentials linear in x and y are locked 
to form the object We will utilize this fact when extracting Fock coefflcients for the Laughlin state in 

the next section. For future convenience we also deflne the short hand notation 

(z) = ( 3 ^ 

_ ^i2Trk^ ^iTTTN^y'^ ^ 
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for the terms entering 0. This allows us to write the basis wave functions 0 as 

Vk (z) = J— fk+N^t (z) . 

V^ tez 


( 4 ) 


In this manner we can isolate the Fourier factors fk (z) that make up the single particle orbitals rjk- Note 
that if we set t = and pk = ^k, then fk (z) can be identified with the single particle eigenstates on a 
cylinder as (k {z) = in the ordinary Landau gauge and physical coordinates. 

2.1. The Fock basis 

In this section we summarize the essential features of the occupation basis, or Fock basis, for many-body 
quantum Hall states. The single particle orbitals defined in equation 0 are labeled by the momentum label 
k. A Fock state is defined as a state with a definite set k of single particle orbitals rjk occupied. We write 
this state as 


5k (z) = A 


Afe 




( 5 ) 


where A is an anti-symmetrization operator over the different coordinates. The anti-symmetrization operator 
is efficiently implemented by constructing the slater determinant of rik^ (zj). A generic many-body state can 
thus be written as 

( 6 ) 


V'MB (z) = Y “kS'k (z) , 


where the sum is constrained to 1 < fci < ^2 < ... < < N^. 

For bosons, the anti-symmetrization operator is replaced by a symmetrization operator S and the 
determinant becomes a permanent, which is harder to evaluate. Also, the set of momentum labels need 
not all be distinct, so the < is replaced by < in the sum over k. 


2.2. The Laughlin wave function 


The original construction of the torus Laughlin state goes back to Haldane and Rezavi |HR 85| by generalizing 
the short distance behavior on the plane to also be true on the torus. The complication on the torus was 
to find a proper center of mass (CoM) piece As {Z,t) that together with the Jastrow factor and Gaussian 
piece would give the desired single particle periodic boundary conditions. The resulting wave function for 
the Laughlin state at = 4 is 


fjs iz) = Af (r) E. ^ I ^ 

Kj A / 


( 7 ) 


where Z = = Zi — Zj. The function "di is related to the generalized z^-function in Q as 


(z| r) = z? 


L 2 J 


[z\ r), and has zeros at z = rz -|- mr. 


The label s enumerates the q different degenerate ground states that exists on the torus |Ha 
states only differ in the CoM function As, given as 


As{Z,T)=d 


n!!-i 


{qZ\ qr). 


These 


( 8 ) 


The q states are all related by rigid magnetic translations of all the particles. As such, the different states 
are transformed into each other as ^y'^'4’s = f’s+i- Note that f’s+q = f’s since under a rigid translation 

of q steps Z —Z + rq^ = Z + t, and z9 is quasi-periodic under this shift. The CoM label s is related 
to the total momentum of the state as Altotai = WsmodA^^, which also shows that there are precisely q 
independent values of s. To specify periodic boundary conditions, s is an (half ) integer when Ng is (even) 
odd. 

The normalization factor Af (t) comes from constructing the Laughlin state using conformal field theory 
and is 

qNe 


A/'(r) = - 


(rY 


71 {tY 


gNe(Ne-l) 


+ 1 ' 


( 9 ) 
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where rj (r) is Dedekind’s 77 -function. This normalization was first introduced by Read |Rea09| and later 
generalized to the full chiral abelian hierarchy in |FHS14| . The normalization coefficient JV (r) has a particular 
T-dependence that ensures that 'tps transforms as 

f’s ^ y^^Ss,s>f’s' 

s' 

s' 

under the modular S transformation r —^ and T transformation r —^ r -f 1 . The matrices Ss,s' = 

24) gj.g |-]^g ixiodular S and T matrices of the CFT for the Laughlin state. 

v? ’ ’ _ 

The modular covariance of ( |10[ ) is related of the redundancy in the parametrization of the torus using 
r. All T related by S and T transformations are identical and this has to be reflected in the transformation 
properties of the wave functions under such transformations. In particular, the space spanned by the g-fold 
degenerate should be invariant. 

The real space wave function ips is however not properly normalized. There exists an extra - unknown 
- normalization factor Afo that can not be obtained by CFT arguments. The normalized Laughlin wave 
function ' 0 z, is related to ips as 

( 11 ) 


and the general consensus is that in the thermodynamic limit (at fixed r) Aq is constant - independent 
of r - and only depends on the particle number Ng. To say that A/q is constant is analogous to claiming 
that the Laughlin plasma |Lau83| is in a screening phase. Several works have added to the consensus of 
screening [RRl 11 IBBR12] but they have not considered the TT-limit. 

From the analysis of the modular properties of we can reassuringly see that the proposed normalization 
(© ensures that the density of ifg is modularly covariant. This means that 


\^S {z,t) \ = \'lps {z,T+l)\ 



which is a also a property that we request from ifr- As a byproduct we note that the unknown normalization 
factor Ao must satisfy the same modular covariance 


|Afo(T)| = |A-o(r + l)| 



( 12 ) 


Although unknown, Ao can be computed numerically for small systems by using that ipr is the 
normalized ground state of the Haldane pseudo-potential |Hal83a| . By comparing the analytical ifs in (|^ 
to the numerically obtained t/jl, Aq can be extracted. In Figure]^ Ao is shown for Ng = 6 particles in the 
entire r-plane. As expected Ao is constant and independent of r in a large section of the r-plane. However, 
when T 2 > 10 then Aq develops clear r-dependence, see Figure This change in Aq signals that we are 
entering the region where the thinness of the torus becomes noticeable. Physically this means that one of 
the torus handles is so thin that the plasma stops screening. We will return to this in Section where we 
will use the analytic Fock expansion developed in Section]^ to analyze the behavior of Ao in the TT-limit. 


3. Normalization and Berry Phases 


We can indirectly gauge the r-dependence of Aq by computing the Berry curvature of Q as a function of 
the geometry parameter r. Computing the Berry curvature is analogous to computing the quantum hall 
viscositv |ASZ951 IRea09| . Just as in Refs. [Rea09| and |FHS14| we make the ansatz that the normalized 
Laughlin wave function can be written on the form 

V’L =A/'o"Vf7/)({z};T), (13) 

where is holomorphic in z and r, and Aq is constant. The power of T 2 is then P = 2 ^. Note that ^ 
has precisely this form. After a short calculation, one finds that the Berry potential Ar = {ipL is 

At = Af = — = — 2^. In terms of and Ar^ this means that 


A't^ — Ay ~\~ Ay —- 

T2 

Ay^ = iAy — iAy = 0 


qNg 

4t2 


(14) 
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Figure 2. Scan of the true normalization J\fo of the Laughlin wave function 0 in the whole r-plane. The 
red dashed lines mark the fundamental domain and the black dotted lines mark the region zoomed in on in 
b). Both panels are for Ne = 6 electrons. 

a) The normalization A/q is independent of ri and depends only weakly on T 2 in a large region of the r-plane. 

b) Focusing on a small region near 0.6 < T 2 < 1.6, small r-dependent finite size modulations of A/o can be 
seen, but these appear only in the 4:th decimal of A/q- These decline at larger system sizes that are now shown 
here. Note how the the pa tter n of modulations conform to the invariance under modular transformations 
T ^ T + 1 and r ^ in ( |l2[ |. 
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b) 


Figure 3. Berry connection in a large part of the r plane, for Ne = 8 particles, for illustrative purposes. The 
black lines show the boundaries of the fundamental domain and a S'-transformed (r^ ~ r) thereof. 

Note that the T 2 axis is logarithmic. 

a) Plot of When Mq is constant one expects a constant T 2 Ar-i = — = — 6 which can be observed 

in a region around r ~ When T 2 ^ oo then AriT 2 diverges. 

h) Plot of T 2 At 2 - Here it is clear that when r 2 > 0 then Ar 2 ~ 0 for most tori. We note however that 
for small values of T 2 then At 2 (for small ri). This region can be understood by noting that it is the 

tS-transformed image of the region r 2 > 5 for • 

Taking both and At 2 into consideration we conclude that the normalization Mo is constant for a region 
around r 5^ 2 but deviates when T 2 ^ oo or when |r| is small. 


We note that As is a gauge dependent quantity and will be sensitive to the r-dependence of the phase 
of 'ijjs- Nevertheless it is a good first test to detect when Mo might not be constant. 

We use importance sampled Monte Carlo integration to evaluate the Berry connection A-^-^ and A^.^ 
numerically as 


At, = lim {{tfs (t) \ips (r + e))) 
e-fO e 

e-^O e 

In Figure we scan Ag over a large part of the r-plane. For illustrative purposes we show a relatively 
small system of = 8 electrons in the figure. The features in this picture are also present in larger systems 
but will occur at different values of t. In the figure we can see that At^ « ~ 4 ^ ^ large region around 

T = i but that it deviates significantly from the expected value when r 2 —>■ 00 . For = 8 particles this 
deviation begins at around T 2 > 5 and is almost independent of ri. 

We also see - Figure]^) ~ that At^ « 0 everywhere except in an approximate circle of radius r « 0.1 
centered at r = zr. This non-zero value can be understood by appealing to modular covariance and noting 
that the this region is the modular image of the region r 2 > ^ « 5 for Ar^ in Figure]^). 

Since the Berry connection is not a gauge invariant quantity, it can be difficult to see what part of the 
deviation in Ag with respect to Ag is due to an actual r-dependence of the normalization, and what is due 
to a simple r-dependence of the phase (at fixed z). For that purpose, we also compute the Berry curvature 
which is a gauge invariant quantity. Under the assumption of a constant Aq the Berry curvature is 

T ——lOT - — f) A —BA _ 

*/ Tir2 — TT — ^T2-^Ti — A 2 ' 

The curvature is closely related to the Hall viscosity |ASZ95I IRea09| which can be used as a probe to 
distinguish different topological phases at the same filling fraction. Hall viscosity is conjectured to be a 
probe containing the same information as the shift 5[Rea09]. On a spherical geometry the Laughlin state is 
characterized by the shift iSgphere = g in the relation 

^ Afg tSgphere 5 
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b) 


Figure 4. a) The viscosity units of the average orbital spin s plotted in the entire r-plane for Ne = 8. Just 

as for the Berry connection A^and Ar 2 we see the expected curvature J^riT 2 ~ ^ large region of 

4 ■’"2 

the r-plane. Note that as TtxT 2 is a gauge invariant quantity the deviation from — when r 2 > r = 3 is 

4''"2 

nicely reproduced in the circle |r — ^ | 

b) The average orbital spin s plotted for r = it 2 and several different Ne- The Inset shows same data as a 
color plot. In the figure it is clearly visible how the plateau at s = | widens monotonically with system size. 
We also see that at large T 2 all system sizes transition to s = ^ which is expected for Fock state. 


between the number of flux quantum and the number of electrons. The shift appears as a consequence 
of the curvature of the sphere and the non-zero orbital spin s of the electrons. When the electrons move 
over the surface of the sphere a Berry phase is accumulated giving rise to a extra effective magnetic flux, 
^sphere = 2.s |Rea08LlRRTT) . 

On the torus there is no curvature, so here 5torus = 0 and = v~^Ne precisely. However, as ^sphere 
is a topological characteristic of the Laughlin state, the same information should exist also on the torus, the 
question is how it manifests itself. Read |Rea09| showed that the topological information could be extracted 
by studying the Hall viscosity (ASZ95] rj^ of the quantum fluid. Read conjectured that the viscosity in the 
thermodynamic limit (on any geometry, but especially on a torus) should be 

= Th^sphere = ^flUS, (15) 

where here n is the electron number density. In a later paper Read & Rezavi [RRll| numerically showed this 
result to hold for the Laughlin state (amongst other things), when r was close to t = t. 

To analytically compute viscosity for a many body state is usually difflcult, but it is simplified 
dramatically if the real space wave function can written as 

V'MB ( 2 ) = tI"/( 2 ; r), (16) 

where / is holomorphic in r. Under this assumption the viscosity can be computed to be rj^ = y |Rea 09]. 

This can be reformulated as the intensive quantity of effective average spin s as 


2P 


(17) 


The Laughlin wave function -(/is in (71 is precisely of the form (161 where Af (r) contains the factor T 2 ^ , 
so s = I for ips- In a similar manner the viscosity for the rest of the chiral Haldane-Halperin |Hal83al IHal83b] 
hierarchy on the torus |FHS14| can also be computed. In doing so a viscosity is obtained that also agrees with 
Read’s conjecture (15l. 

The numerical work by Read & RezayijRRTT] showed that s = | and supports that Read’s conjecture 
holds in the thermodynamic limit for the nearly square tori. This lends indirect evidence to the screening 
properties of the plasma analogy. However, in other numerical studies of the viscosity |ZNS131lFHS14irTH15| 
it is clear that s has significant r-dependence, if the torus is asymmetric enough. 

In Figure |^) we plot the effective shift s = — 2 r-plane for Vg = 8 particles. Again the 

the same features as for At-^ and At^ are present. However, one important difference is that in the TT-limit 
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Figure 5. The Berry phase accumulated by -ipL under a T transform from r to r + 1, as a function of T 2 
for n = 0 and several Ne- Dashed lines show the expected values (j)q- given by ( |18| |. Inset shows the same 
data as the difference <j)q — (fi-j-. The path r (s) is discretized in 100 steps with 7 x lO'^ MC points at each 
step. 

It is clear that agrees well with \is\ when T 2 is not to large. However for large r there is a clear deviation 
between (pq- and (pq- which is not an effect of the discretization of the path. Note that as T 2 ^ oo then 
tpq- —>■ ^ instead of pq- 0. 


s stabilizes at ^ instead of the expected s = |. Note that the features at large T 2 are modularly mapped to 
the ringlike structure at small |t|. 

In Figure Ep) we get further indication that the transition from s=|tos = ^isa generic feature 
independent ofsystem size. In the plot we see that the transition to s = ^ happens for all systems sizes 
examined, but the point of transition scales as T 2 oc y/Ne (see inset). 

This TT-limit value of s = | can be understood if the state at r 2 —t oo is described not by a strongly 
correlated fluid, but rather by a single slater determinant, i.e. a Fock state ifs = 5^k {z) as in (§. We will 
expand on this observation further in Section For now we simply conclude that in a region around r « i 
the Berry curvature has the desired properties. The number of MC data points is here 10® for all system 
sizes. 

Returning to the Berry connection, we can now also compute the Berry phases associated with the 
modular T-transform r —>■ r + 1 and 5-transform r —>■ — y Starting with the T-transform, we find that the 
accumulated Berry phase for a straight path from r to r -|- 1 is 

(fr = [ dri (V'L 19^1 V'i) = = -yr^- (18) 

Jo 4r2 

where we assume that No is constant. 

The accumulated Berry phase (fj- can also be approximated numerically by discretizing the path from 
r to r -f 1 into n steps and computing the cumulative overlap along the path as 


= J ds {if (r) \dsif (r)) « 9 I {ifr {tj) \ipL (Tj+i)) 


i=i 

A comparison of the numerical value (fj- and theoretical value (fj- can be seen in Figure|^for 5 < Ng < 20 
and 1 < T 2 < 20. The numeric data is represented by points and the theoretic expectation is shown by by 
dashed lines of the same colors as the points. One can clearly see that for small T 2 , the numeric phase (fj- 
agrees very well with (18l. It is also possible to see that a t lar ger T 2 the the accumulated phase cfj- saturates 
at (j)j- « 0.5 instead of (fq -0 as expected. In Section 5.3 we will show that this is actually constant is 
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Co 

1-^ 



0.4493 1.0 2.2255 4.953 

^2 


11.0232 


Figure 6. The Berry phase it’s accumulated under an S transform from t to — as a function of T 2 for 
ri = 0.4 and several N^. The dashed lines show the expected values (j)s given by ( |20[ |. The inset shows 
the same data as relative deviation ~ /^S- The path r (s) is discretized in 500 steps with 10® 

MC-samples at each step. 

In the main plot the numeric (j>s follows the theoretically expected 0^ for all T 2 and system sizes considered. 
The inset plot shows no structure in the relative difference between cj)s and 05, and hints that this is 
dominated by numerical noise. The main reason for this noise is that the len gth of the path increases with T 2 
such that the discretization error grows. We note that 05 agrees well with |20| as long as T 2 is in a (large) 
region around T 2 = 1. 


actually f/j- =_^in the limit T 2 —> cx). In the inset in the same figure, one sees that the T 2 at which (fq- 


deviates from (181 grows monotonically with Nf.. This hints that in the thermodynamic limit the transition 
will be at T 2 — f CX). 

We may strengthen the picture that the normalization is constant for t near r « z by also considering 
the Berry phase accumulated under an 5-transform. For the 5-transform one has to be a bit careful, as there 
is no canonical path between r and — ^. Depending on the value of r some paths would be numerically more 
stable that others. We choose the a path that is self-dual, i.e. where t (s) = One such path which is 

particularly nice is 


T (s) = r ( ^ ) = 


(19) 


where r (0) = re*®. Using that Ar^ = — Ar^ = 0 and ^ = —21nr • ti (s) -|- {26 — tt) T 2 (s), the Berry 
phase becomes simply 


4>s = \ ds {-2\nr ■ Ti {s) + {26 - tt) T 2 {s)) At 

Jo 


qNe 

4 


ds —21nr 


Tl ( 5 ) 
T2 (s) 


-b {26 - tt) 


qNe 


(tt - 26) 


( 20 ) 


where we in the last step use that = 0. Again we assume that JVg is constant. 

In Figure we plot the phase 05 accumulated by performing a 5-transform along the self dual path 
(191 starting at r = 0.4 + it 2 . We see that this particular path gives a (fs that agrees well with the expected 


(j>s as long as r 2 < 6 for all system sizes considered. Beyond this aspect ration it is possible to see clear 
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deviations from (20). However, the deviation can be mostly attributed no Monte Carlo noise (upper left in 
inset). That the noise levels are bigger at larger T 2 is only natural as the length path through r-space will 
increase as \t — i\ grows. 

Taking all the data into account (especially from Ag and ^> 1 x 2 ) we conclude that the normalization Aq 
is constant in a large region around t k, i. We also conclude that the region of constant Aq increases with 
Ne- However we also see that if T 2 is sufficiently large, then Aq develops a clear r-dependence and the ansatz 
(13) with P = is changed to P = This new power of T 2 is consistent with being described by 
only one Fock state. 

To gain more understanding of the nature of the transition to this Fock state, we will in the next sections 
expand ips in a Fock basis and analytically send r 2 —> 00 . 


4. Expanding Laughlin in a Fock Basis 

In this section we will rewrite the real space Laughlin wave function 0 in the single particle orbital basis of 
([^, i.e. the Fock basis. We do this to gain analytical control over the coefficients Uk of the Laughlin state. 
Later, in Section]^ we will use this analytic knowledge to make statements regarding the plasma analogy in 
the limit of a thin torus. In Section we will also outline how to perform the Fock expansion for the full 
chiral Haldane-Halperin hierarchv |Hal83al IHal83bl IFrel3| on the torus. 

We begin the next part |4.1| by displaying and commenting on the main results. Readers who are 
interested in the derivation of these results may consult the supplementary material. In the following part 
|4.2[ we comment on the numerical efficiency of the Fock expansion and compare it to other methods of 
obtaining the same information. 

4-.1. Main results and formulas 

In this part we present the main results and formulas. In what follows we will consider only the fermionic 
Laughlin wave function unless otherwise specified. The fermionic and bosonic Laughlin wave functions are 
so closely related that all results here carry over to the bosonic case with only minor modifications. Where 
there are differences these will be pointed out. 

The Fock expansion of any state will in its most general form involve a sum over all partitions of the 
momentum numbers fci,..., with the total momentum ATtotai = cnodN^. For periodic boundary 

conditions ki G Z are integers and Atotai = {Ne — 1) + mNeUiodqNe for m € Z. Since a Fock basis 

expansion could potentially contain all of these partitions it could become a very large sum. A convenient 
way to organize this sum is by changing to a variable offset from ki which we call T,. The relation between 
the two is 


ki =Ti + smodN^, (21) 

where the are chosen such that The variable s is precisely the center of mass momentum 

label in ([^. Because ki is an integer, the index will be an (half) integer when Ne is (even) odd, to match 
s. 

Using the auxiliary variable the Laughlin wave function can be rewritten as 




tps =Af {t) 


Z {T)Y[rJr^+Az^). 


( 22 ) 




2=1 


which makes explicit the single particle orbital basis. Here, pki [zi) are the single particle orbitals from ([^. 
Note that no explicit anti-symmetrization is needed as the coefficients Z (T) are all fully antisymmetric by 
construction. 

As the basis states rjk (z) are labeled modulo the sum over T is finite. The half-integer offset to T 
is a remnant from the expansion of the Jastrow factor. Note that the CoM index s only enters the relation 
between ki and as a constant shift. This shift serves to relabel the basis states of the expansion. The same 
relabeling can be achieved by acting with the ty translation operator as tyTjY. = r^Ti+s, which shows again 
that all the q different Laughlin states are related by rigid magnetic translations. 

The factor Z (T) is anti-symmetric in the arguments T and is related to the Fock coefficient Uk in ([^. 
As it stands (22) is not of the form of (1^. However as Z (T) is anti-symmetric, (22) can be recast on the 
form ^ with the identification Ok = ^/NelAf (r) Z (T). 
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The Fock expansion coefficient Z (T) is given by 


Z(T) = 



qNe 


+ '^jk + 


Tii 


(23) 


i<j<k 


i<j 


where the Ty are for technical reasons (see the supplemental material) anti-symmetric half-integer numbers 


(24) 


Tij = —Tij. The sum over all the is constrained to fulfill 

N, 

TTi = Tij , 

i=i 

such that each particle has the desired momentum T^. 

To obtain the expression for Z (T) one has to expand the 'di-functions in the Jastrow factors in terms 
of Fourier - or momentum - components This makes it possible to extract the total momentum 

of a specific particle Zi as T,j=i Tij multiplying together all the contributions from all the pairs of 

particles. Thus Ty appear naturally in Z (T) as they enumerate all the ways of forming the total momentum 
Ti from all the pairwise contributions Tij. 

In the above formula 


23 


the Tj. 


are offset from integer by | and subject to (24l. 


This is a direct 


consequence of the fact that each -di-function in the Jastrow factor is taken to the power of q. The higher 


power in the Jastrow factors is also noticeable in the factors Z, 


(</) 


The Z''-- 

T, 


(9) 


can be thought of as structure 


factors for the q™ power of a 'd-function. It enters the Fourier expansion of the i?-function as 




(z|r)) = 


(25) 

and was discussed briefly in Ref. |Frel3| . The sum over T is infinite and offset from the integers by aq. 


Details of as well as explicit expressions for it are given in 

ii 


The factor e™Tj in 


23 


Appendix A 


also comes from the Jastrow factor. It is this piece that ensures that Z (T) 
is anti-symmetric with respect to T. The anti-symmetry comes about since an interchange of and Tj 
will only affect the factor e™Tj , Jt does so by sending e^'^Tj ^ Since Tij is a half-integer then 

Tij — Tji S 2Z -I- 5 is an odd (even) integer and g-^Tj = (—1)“^ 

See Figure 7 d for and interpretation of the sum {^Tij + Tjk + Tk^ as sum over triangles with 

corners in and 

I Ne 

( 2iV TT^ \ ^ 

—— ) that 


come from 


the CFT normalization and normalization of the r]k respectively. 


4-.2. Numerically evaluating the Fock coefficients 

In the previous part we showed that the Laughlin state Q could be expanded in a Fock basis with coefficients 
Z (T) given by (23). Unfortunately the given expression is not particularly helpful when it comes to numerical 


evaluation. The culprit is the simultaneous infinite sums over all of the T^cs. To alleviate this problem a bit, 
Z (T) can be rewritten such that the sums can be performed in a recursive manner. This will substantially 
reduce the scaling of the computation. The details of this manipulation can be found in the supplemental 
material. 

A word of caution should be given. Although the optimized expression for Z (T) can be put on a 
computer, the efficiency in evaluating it is still inferior to that of simply diagonalizing the Haldane pseudo¬ 
potential Hamiltonian. Using exact diagonalization the coefficients for around 12 particles can be extracted 
numerically, but in our current implementation of Z (T), only 6 or so particles can be achieved. 

Our method can as mentioned above not compete with exact diagonalization for the Laughlin state, but 
that is on the other hand not the main purpose of the expansion. An obvious advantage of our approach is 
the explicit access to the r-dependence of the Fock coefficients, something the exact diagonalization can not 
provide. 
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Figure 7. a) A tabular view of relationship between Tij and T; in equation ( |24[ | we well as the balance 
condition Ti = 0. b) How the signs of the terms in Tij are added to form l |23| . In the figure it is assumed 
that i < j < k. 


Also, the reader should also be aware that the pseudo-potential trick only exists for the Laughlin state. 
For states higher in the hierarchy there is no local Hamiltonian for which these state are the exact zero energy 
eigenstates. In these cases the results in this article is the only way we are aware of to analytically extract 
Fock coefficients. 

There is also another setting in which analytical knowledge of these coefficients may be of use. 
This is when considering the non-chiral extensions of the chiral HH-hierarchy. For this extension the 
iF-matrix |WZ91j is spit into two separate pieces as AT = k — R. This can be done even for the Laughlin 
case by splitting K = q as k = q + p and R = p. On the torus, the CFT methods introduced in Ref. 
|HSB~*~08| and Ref. |FHS14| can be generalized to handle general AT-matrices in an analogous way to the 
procedure on the Dlane |HCJV07bl IHCJVOTa] . The wave function that is obtained |Frel3| contains factors 
(ziji {zij\f)^ such that it is not only in the LLL any more. This wave function can be projected 

on the LLL analytically using a basis of coherent states [EVel3| . In this projection, the coefficients Z (r) and 
Z (f) of the underlying Laughlin states appear naturally as ingredients. The details of the analytic projection 
mentioned in the above paragraph will be the subject of a future paper. 


5. Asymptotic Behavior in the TT-limit 


In this section we will return to the main equations (22 1 and (231. We will analyze these expressions in the 
limit of r —>■ zoo (or A —)■ 0 for constant area), which is what we call the thin torus (TT) limit. We are 
interested in studying the TT-limit to answer fundamental questions regarding the plasma analogy. Using 
techniques from conformal field theory [Rea09| a candidate for the normalization Af (r) of the Laughlin state 
can be obtained - see equation ([^. The assumption that Af (t) captures the full r-dependence is related to 
the plasma analogy introduced by Laughlin |Lau83| . It states that the normalization of the Laughlin state 
is the Boltzmann weight for a single component plasma of charged particles with a logarithmic electrostatic 
interaction. As this plasma is know to be screening, the braiding properties of the quasi-particles can be 
deduced, assuming they are separated far enough. On the torus, the assumption of screening properties with 
respect to r allows us to compute e.g. modular group representations. In this section we will check if the 
screening holds also in the TT-limit. We will expand on this discussion in part |5.2| 

we analytically take the TT-limit and extract the 


This section has three parts. In the first part, 5.1 


dominant Fock contributions. In this limit only a few of the Fock coefficients will remain nonzero, and their 
analytical expressions are drastically simplified. We also extract their relative scaling as a function of r. 
This knowledge will then be used in parts and |5.4| to make statements about the corrections to the CFT 
normalization and the Hall viscosity in the TT-limit. 

Section 5.2 concerns the true normalization of Q in the TT-limit. We will use the fact that there is only 
one particular configuration with non-zero Fock coefficient in the TT-limit - the root partition - to compute 
the full normalization when r -A loo. From this knowledge we compute the asymptotic correction to the 
normalization Af (r) proposed by the CFT construction in Ref. |Rea09| . The result shows that the plasma 
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analogy does not hold in the TT-limit. 

Under the assumption that Af (r) is the full normalization, i.e. that the unknown piece Ao cont ains no 
T-dependent contributions, the Hall viscosity r]^ for ifr can be analytically extracted. In part 5.4 we will 


analyze how the viscosity changes as we approach the TT-limit. We will argue that since it is only the root 
configuration that is present in this limit, the viscosity is trivially r]^ = ^hh instead of = ^hn which is 
the value for the square torus. This is the same viscosity as that of a single slater determinant state. To make 
the claim more robust, and also give some indications of why the viscosity is independent of system size in 
the regime of small L, we compute the viscosity analytically for the root configuration and the sub-dominant 
Fock contribution in this limit. 


5.1. Asymptotic scaling of the Laughlin state in TT-limit 

Let us now take a look at the TT-limit for the Laughlin state. In this section we will show that the root 
configuration is the dominant configuration lSSSTl IBKO^ . and we will also give the relative scaling of all the 
two particle squeezed states. 

We start from ( [2^ and use the short hand notation A (W) for the Gaussian weight in Z (T): 

H(W)= ^ (26) 

i<j<k 

We are interested in the behavior deep in the TT-regime, when r —^ zoo. In this regime the Fock expansion 
will be exponentially dominated by the configuration of Tij that can maximize the weight 

(Ty) 


W 


(27) 


i<j 


For large r it is straight forward to show that (z| r)'^ « 1 for — ^ < 3 (z) < ^ and as a consequence 

1 (and zero otherwise) when Tij = ±|. 




2 in this limit. That is, has its largest value Z 

’ 2 J if T, 


(A 


We take as an ansatz that this configuration will also give the maximum for W (^f'ij'^. Explicitly, we choose 
the ansatz 


= I - 0.,,) , 


(28) 


where 9 is the Heaviside function 

1 i< j 

i > j 


= 


In words: T^ = | 


0 

when i 


> j and Tij = — | when when i < j. This particular ansatz is favorable 


y ~ 2 

since it produces the correct root partition. The electron with label i will (when s = 0) have momentum 
ki = Ti = J2j Tij such that ki = qi — ^ (W -I- 1). The value of A (W) for this configuration is 

e (|+i-|)'= i: ^ = |^(iv=-3iv.^+2iv.). 


i<j<k 


i<j<k 


The sum = J2i<j<k ^ = g (-^e ~ -I- 2W) is computed either as the third order polynomial that has 


roots Pq = Pi = P 2 = 0 and P 3 = 1 or by using Faulhaber’s formula ~ \ -|- n). 

The reference weight is thus 

Wtt = (02.2 - 

yields a global minimum (up to permutations of the electrons) 


(29) 

We must still ensure that the ansatz ^ 
of (271, but here we only check that it is locajf] To do so we investigate changes to A{Ne) as one of the 
Tij is varied. Due to the permutation symmetry in A^N^) we can without loss of generality choose i < j. 


There is a subtlety here; as we have chosen the reference TT-configuration to be that of (28 1 , the ordering 


of the electrons is implicit from 1 to W- Another choice of the signs in (28 1 would lead to a corresponding 


permutation of the particles. Mathematically this shows up in (26 1 in such a way that Ty is negative in some 
terms and positive in some others. We can interpret a change Tij —> Tij -\-m as electron i moving m steps to 
the left and electron j moves m steps to the right. Thus - because of the implicit ordering - (for i < j) for 
TO > 0 is a squeeze and to < 0 is an anti-squeeze. 


t A global proof has not been attempted. 
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We start by considering the simplest case of changing only T 12 T 12 +m. The effect on W [Tijj comes 

from two factors. The first is the change in A that is AA = (TVe — 2)m(m — q) and the second is from . 
This factor becomes 



(30) 


for I < m < |. 

By looking just at AA we see that A^l < 0 when 0 < m < q, which translates into the particles 
being squeezed. This must be a pathological result as it would render the TT-state unstable to clustering of 
electrons. However, when we also take into account we get 

Tij 

^ ( |m| ) (31) 

The combined difference AW is 

Air = -2m» + 2W,xR__, , (32) 


which for \m\ < f is always a positive. Note the somewhat counterintuitive behavior that ... OlotoOOOlolO ... 
is more penalized than ... OlOOOttoOOlO ... in the TT-limit. The former would have an excitation energy 
of /S.E = 2 A 3 + A 4 + • • • and the latter Aif = A 2 + 2 A 3 + • • • which is larger. The terms A^ are greater 
than zero and can be though of as pseudo-potentials, although their origin is slightly different. For simple 
potentials like the Coulomb potential Ai > A 2 > A 3 > ... > 0. See e.g. equation (19) in Ref. |BK08| 
for details. Of course, using perturbation theory, the first correction would not be given by the energy, but 
by the amplitudes of the hopping terms. In this case the hopping term for squeezing V^^-i is larger that 
anti-squeezing V 34 . Thus, the relative scaling of the squeezed m = -1-1 state is qe ™'^whereas it is 


anti-squeezed state. 

In Figurej^we see the relative scaling of the root partition in comparison to the sub-dominant Fock-states 
for TVe = 6 particles. We see that the two leading sub-dominant states are a squeezed (OOltoOOlOOlOOlOOlO) 
and a doubly squeezed (OOltoOOlOOOttoOlO) state. The order and scaling of these is well captured by the 
leading order expansion ( [3^ . In the Figure it is possible also possible to identify a triply squeezed state 
(OOltoOOOltoOOintoO) and this is also well descried by (321 taken three times. 

Observe that as the area is fixed to LfT 2 = 2TrNf,q we can rewrite the scaling difference as = 

e where the scale depends only on L, and ri only changes the relative phase. This shows 

that the relative scaling in the TT-limit is independent of system size and depends only on L instead of 
r. This is actually to be expected since in the thin torus limit the only relevant length scale is the short 
circumference L. 

We can now also consider the more generic case where particle i and particle i + A are squeezed m steps. 
The change in A this time is 


AA = (Ne — 2) — mq (W — 2A), 


where A = j — i. Taking into account the change from Z^"^ gives again (|31| but this time with 

J- in 


AW = — 2mf -I- 2 |m| q x 


A TO > 0 

Ne — A TO < 0 


(33) 


(34) 


Putting A = 1 yields the result in (321. Just as with the special case T 12 , AW > 0 for small values of to. 
Note also that the formula for to > 0 can be mapped on the case to < 0 by letting to —)■ —to and A —)■ Ne — A. 
This symmetry is natural since an inwards squeeze of to steps a distance A away can alternatively be seen 
as outward squeeze of two particles Ng — A steps away. 

Looking again at Figure we see that we can identify the one particle squeezed states 
(OOioiotoOOlOOlOOlO) corresponding to (A = 2, to = 1). We can also see (OOlOlotoOOlOOOltoO) which 
is a combination of (A = 2, to = 1) and (A = 1,to = 1) acting on different pairs. The scaling of both of 
these are well described by (341. 
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Relative strength of Fock coefficients In TT-limit 
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Figure 8. The relative scaling of the difference Fock coefficients for Ne = 6 particles compared to the root 
TT- configuration ...010010010010.... The largest coefficients are labeled with the pattern of occupied 
orbitals for that state. The black dashed lines correspond to analytic scaling given by l |32[ | and | |34[ |. The 
analytic expansion fits well with the numeric data. The red text highlights the two orbitals an electron has 
moved between. These states are obtained by squeezing electrons that are 5 =3, 3 + 3, 6, 3 + 3 + 3, 3 + 6 
orbitals away. Here 5 is the number of orbitals between the squeezed electrons and all squeezes are one step 
inward. The notation a + 6 (a + 6 + c) means that two (three) different pairs distances a and b (c) orbitals 
away are squeezed. The analytic scaling is obtained by treating the two (three) squeezes as independent. 


5.2. A/q in the TT-limit 


In this part we will investigate the possibility of TT-limit corrections to the normalization of the Laughlin 
state proposed by the CFT construction in Ref. |Rea09] . If TT-corrections turn out to be present in the limit 
Ne —>■ oo, it will have immediate consequences for the validity of the plasma analogy. The plasma analogy 
states - in the form relevant here - that the free energy of a single component plasma in two dimensions on a 
torus does not depend on the geometry of the torus, provided the area is held constant. The original plasma 
analogy formulated by Laughlin states that the free energy does not depend on the positions of test charges 
in the plasma, provided they are sufficiently far apart |Lau83j . In both cases the validity of the analogy relies 
on the property that the plasma should be screening. 


At the end of Section (2.21 we discussed the behavior of the true normalization A/q for a small system 


of A^e = 6 electrons. We found by consulting Figure that as long as the torus is reasonably thick in both 
directions then A/q is also constant. For r ^ i, variations of A/q appeared only in the fourth decimal (see 
Figure If)- However, looking at more asymmetric tori (t 2 > 10) we noted that Ao deviates from being 
constant and develops clear r-dependence. This we interpret as one of the torus handles being so small that 
the plasma no longer can screen in that direction. As a consequence the plasma stops to screen and the 
normalization obtained from CFT is incorrect. 

We will now develop an approximate expression for the r-dependence of Aq in the TT-limit. This is 
enabled since the expression for ijjL simplifies in this limit. Due to the differences in scaling, only one of the 
Fock states will survive and 'ijjL will be just a product state 


V’L —^ ^TT {z) — IJtT (z) ■ 

This product state is precisely that of the root partition. Also, since ipL = A+^V's 
formulated as 


(35) 

the limit can be 


V'L ^ AAq-W (r) Z (T) V+el^TT {z) , 


(36) 


where both Af (r) and Z (T) are known analytically. Comparing (351 and (361 leads to the the form of Aq in 
the TT-limit (t 2 —> oo) as 


Mo Mtt =M{t)Z (T) v+e!- 


(37) 
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Laughlin normalization 



Figure 9. Normalization of the Laughlin u = ^ wave function defined in |7| for rectangular tori, a) 
Com parison of the true normalization A/o (solid lines) to the analytic approximmion A/tt (dotted lines) in 
( |39| b) The quotient A/q/A/tt as a function of T 2 . When T 2 > (marked by ’X’) then A/tt ~ A/q to very 
good accuracy. 


Note that for q = 1, A/q = Att for all t as in this case there exist only one Fock state, the filled Landau 
level. 

The expressions both for Af (r) and Z (T) do also simplify when r —>■ loo as all sub-dominant terms 
will vanish. For instance, the Dedekind eta function rj (r) that appears in Af (r) can be expressed as 
77 (r) = (1 — so the a asymptotic expansion is simply ry (t) —>■ Inserting this 

into ([^ gives asymptotic behavior of Af (r) as 


Af (t) 


oiTrr i [(3iVe(Are-3)+2] ' 


Similarly for Z (T) the dominant r-dependent factor of (231 is given by (291 and so 


Z(T) 


2fV^ 

T2 


^^^r^{N^-3N,+2) 


in the TT-limit. Inserting the above simplifications into (371 gives 


A/tt = \//^ {‘^qNeTT^) T2 


(38) 


(39) 


We should at this points ask ourselves how well Att actually approximates Aq, as it is obtained from 
a rather crude analysis. It turns out that A/tt does a good job, provided r is large enough. In Figure is 
displayed a comparison of Aq and A/tt from r 2 = 1 and deep in the TT-regime. The solid lines display Aq 
and the dotted lines A/tt- Two things should be noted in this figure: 

Firstly, from Figure it was already clear that while Aq is constant over a large region of r-space, it 
does depend on r for thin tori. In Figure]^ we now see that when this happens, e.g. at r 2 « 6 for Ng = 7, 
A/tt matches Ao very well for all values of Ne plotted. We thus conclude that A/tt manages to capture the 























Success and failure of the plasma analogy 


for Laughlin states on a torus 


19 


a) 


b) 


Quotient of Nq and Njj 



Scaled Q as function of L 



Figure 10. Comparison of the modified TT-normalization A/iJirp ( |40| and the exact normalization A/q of 

the Laughlin v = | wave function at rectangular tori, a) Plot of Q = Afo/Aff.r^ as a function of T 2 . The 

1 

revised normalization is accurate to within factors of size unity, b) The scaled quotient Q as function of 
the shorter length scale L. This shows that A/o and A/ij.rj, only differ (up to small finite size deviations) by a 

size independent factor that only depends on the shorter length L as (L). The dashed line (L = tt-\J^ q) 
shows the break point between the two regions of A/’.^rp. 


full T-dependence of Ao where this exists. We should however note that in the region where Ao actually is 
constant, then the approximation A/tt breaks down as it is not constant there. 

Secondly, the maximum of A/tt (at T 2 = lies in a region where where Ao is still approximately 
constant. Under the assumption that this will be the case also for larger system sizes, we may refine the 
approximation of A/tt in the region 1 < r 2 < . We simply assume that if T 2 < then no r-dependent 

correction is needed and if T 2 > then the correction should be according to (391. The constant correction 

^ 37Vc 


is naturally set to the value of 

A/tt Ei-t 

T2 

_ 3jVe 

TT 

[ A/'tt(t = *^) 

1 < r2 

< 

Me 

TT 

A/pT = N 




[ A/pT {t) 

T2 

> 

Me 

TT 


r 

N,- 

5-1 

4 


V CTT / 



1 

i 

-1) 

0—1 

gZTrrL- 


(40) 


1 < r2 < 


T2 > 

^ - -TT 


We define Q = -^P- as the quotient between the exact normalization and the revised approximation Aq-p. In 
Figure 10 1 we show that the revised normalization App is correct up to factors of order unity for all t. To 


the best of our knowledge this is the first estimate of the normalization of the Laughlin state spanning both 
the thick torus and thin torus regimes, although we acknowledge that a calculation for the thick cylinder has 
been performed in Ref. |TGV09] . We also note, in the same figure, that Q is close to 1, but often slightly 
bellow, for all system sizes considered, e.g. 1 > Q > 0.8 for Nf. = 7. 

Further, as we argue that L and not t is the natural variable to describe the physics, we may re-plot 

. 4.44^1 


Q as a function of L. This is a reasonable action since T 2 = is equivalent to L = tt 

which is independent of N^- This means that when L < then Ao has r-dependence irrespective of the 

system size. This is another example showing that the physics in the TT-limit is dominated by the short 
torus length, and insensitive to the overall size of the system. 

Figure 10 d plots as a function of L, where we take the root to cancel the difference in scaling 
’—’ _x_ 

at different Ne- As should be apparent, all the curves [Q {Ne, L)] have collapsed approximately on a single 
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Figure 11. Monte-Carlo estimate of normalization of the Laughlin u = ^ wave function defined in Q at 
T = i, for Nf, ■ 10^ configurations. (Inset) Comparison of t he MC normalization A/mc (red) compared to the 

exact A/o (blue) and the approximate (green) in |4o[ |. (Main plot) Plot of / (Afg) = Q estimated by 

MC (red) compared to the exact (blue). This data suggests that / (L) approaches a constant when L —> oo. 
Note that as the MC calculation has a skew distribution, the more accurate (lower) values have a larger 
error. 


curve / (L). Here, / (L) is an unknown function that is independent of Ng. By graphical inspection the form 
of / (L) seems to converge for system sizes as small as Ne = 5. Also the value of / (L) seem to stabilize 
around L « 8£b to / (L) « 0.973 (3). 

To obtain the numerical values for A/q we used exact diagonalization of the Haldane pseudo-potential 
to access the properly normalized ipr. This method will only work for system sizes where we actually can 
perform the exact diagonalization and generate real space samples. To see that Q {Ng, L) « (L) beyond 

this regime we resort to estimating the norm by Monte-Carlo (MC) and importance sampling. We then 
generate samples according to the probability distribution p (z) that we choose to be the absolute square of 
the ly = 1 wave function. The MC estimate of the norm is then 

^MC 


■^MC ~ {‘2'TrqNe 


\Ne 


1 

Z 


E 

i=i 


\^8 {Zj)\ 

p{zj) 


where Z = and Nmg is the number of MC samples. The factor {27rqNg)^‘ is gives the volume 

integrated over, as importance sampling will only compute the average value of 

we compare A/mc to s-i^d Ay^Ti i^or t = i (square torus) and up to Ng = 20 electrons. 


In Figure 


11 


Exact numeric values for A/q exist for up to Ng = 11 but beyond that it is hard to evaluate in real space 


the ground state wave function ipr. In the Figure 11 we see that A/mc and A/q agree (as they should) 


to within MC errors for all system sizes where we can compute A/q exactly. In the main figure we see that 
Q {Ng, L) « f^^ {L) seems to hold beyond where we have exact data for A/q. We note that A/mc has a tendency 
to overestimate A/q, which is why the lower values {e.g. Ng = 19) has a much larger error than the higher 
ones {e.g. Ng = 18). We conclude that using Monte Carlo strengthens the picture that Q {Ng,L) « (L), 

and that f {L > 8£b) ~ 0.973 is constant. Another way of stating the same is that at r « z is correct 
up to a factor 0.973^', and sub-leading multiplicative corrections must be polynomial in Ng. 


5.3. Modular Transformation Phases in the TT-limit 

We are now in a position to answer what happened with the Berry phase (fx in the thermodynamic limit 
T 2 —t oo in Figure In the TT-limi t we kno w th at ipr is just a slater determinant, so we know that it is 


on the form used to obtain equation (|18|) and (|20|, but his time with P — However we also know that 


tjjs needs to be normalized by the factor A/tt iu (39) which contains a factor of 12 . As a result, the 
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Figure 12. plotted as function of L, such that all data collapses on the curves — (dotted) and 

T ^ i ^ 

— =7 -? (dashed). The vertical (dotted) line marks L = n-J -^q . 


the normalized tpf = 
state is therefore 

ATT) _ _ 


— —A— actually has the form rf e ’12 [{z} , t). The Berry phase picked up by this 


Ne q-1 

'rri - TT . 

^ 4t2 12 


(41) 


Now, using the relation L^T 2 = 2'KqNe we can rewrite the two curves (181 and (411 as 

,rr.rr., A 

and 




247r 




q-l 
~ ^ 12 


Thus plotting (pT in units of L instead of T 2 all the data for different Nf^ should now collapse onto the same 
curve. Indeed, looking at Figure 12 this is precisely what we see. Also, the crossover from to 

happens again around L = 7rw|g. 


5.4- Hall viscosity in the TT-limit 


In this section we will investigate the consequence for the Hall viscositv |ASZ95| when there is r-dependence 
in A/q- From studying px in Figure (12) we know that the proper parameter that determines pT is L instead 
of T. It would thus seem likely that L would also determine the behavior of the viscosity s. 


An example of this is shown in Figure 13 where the data from Figure (|^ is replotted as a function of 
0 < L < y/2TTqNf,. The different curves extend to different L, since the torus area increases with system size. 
Two things should be noted here. Firstly, all the curves of s for different are approximately identical 
when plotted against L. This again lends support to the interpretation that it is L that is important for the 
physics. Secondly, as the thin torus is approached - L —>■ 0 - the value of s drops to s = 4 from the square 
torus value of s = |. This indicates that there are extra Berry phase contributions from the normalization Aq 
and is direct evidence that the plasma is not screening in the TT-limit. This last observation leads naturally 
to the need for systematic corrections to Ao, as a function of r. 


Using the results from section 5.1 and equation (40) we can approximately compute the viscosity of 
the Laughlin state in the thin torus region and in the square torus region separately. We first consider 


V' = 


where A/’^-p is constant for L > 7r-\ 


|( 7 , giving s = I as above. However for L < 


\q there is 


an now an extra contribution r. 


^( 9 - 1 ) 


Thus, according to 0 the viscosity will not be s = | as expected 
from the plasma analogy. Instead we see that the viscosity is s = | — = 1. This is the same viscosity 

as that of the fully filled LLL, or more generally any product state of LL orbitals. In hindsight this result is 
simple to understand physically by remembering that the Laughlin state is the exact zero energy eigenstate 
of the Haldane pseudo-potential Hamiltonian |Hal83a| . In the extreme TT-limit all hopping elements vanish 
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Ne=3 

Ne=4 

Ne=5 

Ne=6 

Ne=7 

Ne=8 

Ne=9 

Ne=10 

Ne=ll 

Ne=12 

Ne=13 

Ne=14 

Ne=15 

Ne=16 




• Ne=3 

• Ne=4 

• Ne=5 
o Ne=6 

• Ne=7 
» Ne=8 

• Ne=9 

• Ne=10 

• Ne=ll 

• Ne=12 
o Ne=13 

• Ne=14 
» Ne=15 

• Ne=16 


Figure 13. (Data Points) The viscosity (measured as effective average spin s) of the Laughlin state as 
function of L for different N^. Note that the curves are almost identical for all Ne even for as small systems 
as Ne = 3, hinting that the physics is dictated by the shorter length scale L. For L > 10£b the viscosity has 
stabilized at s = |, which is the expecte d th ermodynamic value. 

(Solid Lines) Approximation of s using \iA\ for different values of Ne as a function of L. Note that this 
function has the right features of rising from s = | along the same curve almost independent of Ne- For 
L < 5 the analytic and numerical results in the plot are in good agreement. For L > 5 the curves start to 
deviate as the approximation breaks down. 


and the Hamiltonian is purely electrostatic. The Hamiltonian is therefore diagonal and the ground state has 
only one non-zero Fock coefficient; that of the root partition. This wave function (z) has precisely the 

Ne 

non-holomorphic factor T 2 which gives s = | ■ 

For completeness we mention the case of g = 1, which is an unnormalized ly = 1 wave function. In this 
case A/q = A/tt exactly as there is only on Fock state. Thus 


Afo=Af(r)Z(T)\/Ag 


(2A^) 


fVe(-^e-3) 

77 (r) 


+1 


{TijGZ-l-i} 


which is purely holomorphic as all T 2 factors in Z (T) and A/'(r) cancel. Thus we see that g = 1 has s = ^ 
for all values of t - as it should, since it is only a single slater determinant. This is also an example that 
shows that to compute the Hall viscosity analytically it is essential to have control over the non-holomorphie 
T-dependence of the wave function. None of the holomorphic factors will affect the value of the viscosity, 
which allows for the possibility that A/q can have holomorphic dependence on r without affecting the viscosity 
of the state. 


5.4-1- Approximating the TT-limit viscosity In section [5.1| we computed the relative scaling of the root 
partition and the one-pair squeezed states. In this section we will now compute the viscosity for the Laughlin 
state in the TT-limit under the approximation that only the root configuration and the leading one-pair 
squeezed states contribute. We will obtain a viscosity that is not dependent on ti, and where L instead of 
T 2 is the relevant length scale. As the approximation is quite crude, we will not be able to reproduce the 
full features of SeS with this analysis, but we will see good agreement when L < 5. Due to the difference in 
scaling between the root partition and the sub-dominant contributions the viscosity will reduce to Seff = 5 
in the TT-limit. 
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We thus consider the root partition of the q — 3 Laughlin state \TT) = ... 001001001... as well as the 
sub-dominant configuration + = .. . OOlOOOttoOOlOOl .... The number of nearest neighbor one-pair 

I 4 

squeezed states are W- From (321 we get the relative scaling between these terms as a (t) = . The 

normalized asymptotic state can thus be written as 

|r) = N (r) I^ITT) + « M E I*’ * + (42) 

where N~'^ (r) = 1-1- Ne9e~^^^^ is the normalization. To compute the viscosity we use 
1 2r^ 

S — — % ^ ^ ) (^Ti^k ■ ^r2^k ’ '^ri^k) 5 


which is an (for s instead of rju) adapted version of equation (2.99) in Ref. |RRTT] . When applied to 
|t) and identifying qq = N (r) and Oi = N (r) a (r) we obtain 


-1 1 ' 


N, 


ipr^N ■ dr^N — dr^N ■ D^N) 


+ (Na*) ■ dr, (Na) - dr, {NaP ■ dr, (Na)) 


i=l 


— z2r| {dr 


(Na*) dr, {Na) - dr, (Na*) ■ dr, {Na)) 


8'7rT2 

64e Tr’^rf 


23047r4e^ 


2 ’ 


(43) 


(44) 


(e^ + W) ^ (e^ + Ow) 

which apart from a factor of Ng in the denominator only depends on L. The contribution from the first term 


(oo) in (42 1 disappears completely as dr,N = 0. Note that s (t) does not depend on ti and only on L. This 
is an encouraging result as it tells us that the viscosity is insensitive to skewing of the torus geometry. Again 
a natural result if we expect the physics to be dictated by the shorter length scale L. 

In Figure 


13 


independently of W 
which starts to deviate from s = 


we plot (44) as a function of L and N^. We observe that s (r) rises from s = ^ almost 
If we set (artificially) W = 0 we get the idealized upper curve s ( t ) = ^ -|- 230471 - ^ rt- 




at about L > A. We can with this construction show explicitly that in 


0 when L —)■ 0. 


the deep TT-limit we will approach s = ^ as p 

Now, if we compare to the exact viscosity in the same figure we see that there is good quantitative 
agreement between (44) and the viscosity of the full Laughlin state when L < 5. For L larger than this 


value the analytic approximation for the viscosity breaks down. We stress that we should not expect perfect 
agreement for all L with this simple model. The reasons are several. Firstly, when r 2 —?> 1 there are going to 


be increasingly more states that are relevant, which we ignore here. Secondly, the relative scaling in (34) is 
only the leading expansion of relative scaling, and will also break down when r 2 —>■ 1. 

We summarize this section by stating that the viscosity in the TT-limit is not s = | as expected from 
the Read conjecture. However, the L-dependence of s is - up to small corrections - independent of system 
size. Numerical studies of e.g. v = suggests that this holds true for, at least some, more complicated filling 
fractions |Morl5| . 


6. Fock Expansion of the Hierarchy States 

In this section we will generalize the results obtained in Section for the chiral Haldane-Halperin hierarchy 
states recently constructed on the torus using CFT techniaues |HSB~*~08l IFHS14| . The Laughlin wave 
function fits in to a much larger framework of trial wave functions, namely the Haldane-Halperin hierarchy 
of incompressible states |Hal83al IHal83b| . These states have been realized using correlation functions 
in CFT on both plane jHCJVOTbl IHC.TVOTa) . sphere |Kvol3| and torus [HSB~*~08[ IFHS14| . In order to 
make this section reasonably self contained we review the basic structure of the torus Haldane-Halperin 
hierarchv |Hal83al IHal83b| wave functions. We here summarize the construction of the Haldane-Halperin 
hierarchy on the torus and refer to Ref. |FHS14| and references therein for more details. 
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The topological information regarding the quantum Hall fluids are encoded in the Wen-Zee 
it'-matrix |WZ9l] . The matrix is symmetric with integer entries and has the same dimensionality n as the 
level of the hierarchy state. The itT-matrix formalism describes n different layers with particles in layer a. 
The distribution of particles in the different layers are determined by the requirement that the quantum Hall 
droplet should be homogeneous. Under this requirement, the group sizes are given by = N^- 

The Ailing fraction is readily extracted as u — ^a /3 |t| 

For the Jain series pai07| at v = , the /f-matfix is u-dimensional and has a particularly simple 

form Kajs = m + Saj 3 - The Jain series is a subset of the chiral hierarchy, which will be considering in this 
work. As examples, the AT-matrices for the Laughlin v = - states and the v = i Jain state are given by 


K = q and K = 


3 2 

q and a = [ ^ 3 

Electronic wave function for the HH-hierarchy are constructed using correlation functions in CFT. The 
correlation functions consist of insertions of n different types of electronic operators, one for each layer, and 
a neutralizing background charge operator. Many-particle wave functions with well defined single particle 
boundary conditions are identified with the conformal blocks of primary operators in the CFT. At filling 
the many-particle wave functions can be written as 


fraction v = ^ 
q 


(^) = Af (r) li ^ I """ {Z, r). 


(45) 


i<j 


should be interpreted as Kap for the groups a, f3 that i,j belong to. Just as with the Laughlin 


Here AT, 

state, the label s enumerates the q different degenerate states that has to exist when the denominator of 
V = - is q. Also, in analogy to the Laughlin case, these states can be transformed into each other by rigid 
translations of all the particles. The CFT construction also gives a proposal for the normalization factor 
N{t). 


Note that for two particles in the same group, the anti-symmetry properties of (451 are dictated by the 
Jastrow factor ''^1 (zijlr) alone and is not influenced by the CoM piece (z,t). This means that 
each group is fully anti-symmetric within the group. The same is not true for particles in different groups. 

The chiral hierarchy has the property that all the edge modes travel in the same direction. This manifests 
itself in that all the eigenvalues of AT are positive. This enables us to parametrize AT using charge vectors qo,. 
These are real n-dimensional vectors such that Qq -q^ = Kap. The vectors q^ span an infinite n-dimensional 
charge lattice F. The charge lattice F is important for the construction of the CoM function. The CoM 
function is given by 

(46) 


{z,t) — g“'^(q+hd^g*27r(q+ha)(Z+t,)^ 


qer 


where Z = j^^^fliZi. The parameters and tg carry the many-body momentum index. They are 


chosen such that for periodic boundary conditions e 

„z27rs^ / ^\p{N(j,-\-Kii) 


„227rhs 


= (- 1 ) 
^N, 




z27rhs 


and e ”0 = 

e 1 In the last equation appears the total charge Q = ' qi of the combined charge 

vectors. To ensure that all the particles are at the same flux the homogeneity condition 


Q 


N, 


i=i 


= N^ 


(47) 


has to apply. This can also be formulated in terms of K^p as the homogeneity condition that dictates 
the group sizes; J^p KapNp = Np. To obtain the physical wave function from the CFT conformal blocks, 
equation (451 has to be anti-symmetrized and external translation operators have to be added as 


V's=-4 < n 


l)“-i 




(48) 


,a=l 


where the A stands for anti-symmetrization. The operator is the torus counterpart of the external 
derivatives on the plane gj anti-symmetrization in the above equation is there to 

ensure that all particles are indistinguishable. For the Laughlin state in ([^ the A is redundant as that 
state is by construction an anti-symmetric single-component wave function. Note also that A acts primarily 
between groups as all the particles within the same group are by construction anti-symmetrized. For a 
detailed discussion about the construction of D(q), ifg and ifs we refer to |FHS14| and references therein. 


t In this article we are exclusively working in the basis where t = (1,..., 1). 
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The Fock expansion for the chiral hierarchy in (451 can be obtained by straight forward generalization of the 


method used to analyze the Laughlin state. As such, also the form and interpretation of the different terms 


are very similar. Thus, the more general (451 can be written as 

Afe 

f’s=^f{T) 

{T,ieZ+^} 


-^Ve 


(49) 


The structure is again analogous to that of (22). The main difference is that there is an extra summation 
over the n indexes rua, a = l,...,n. These indexes TOq, come from the CoM function (461 and are here 
brought to the form where the they only appear in the definition of the momentum ki 


ki — T Nfp ^ s. 


(50) 


As (49) is written, the sums over ma can immediately be performed to produce X){m } IIq; CfcjVa = 

(zjv„). However we write (491 with this summation unperformed as it only affects n of 


-) " UcVk^ 


Ti + Np which means that we can not write (49) on 


the Ne factors in HE Cfei (zi)- 

This time, Z is not invariant under 

the form X){Ti} 11^ fi Iki {zi) with a simple expression for Z (T^). This is not to say that it cannot be 

done, only that Z (T^) will not have the simple structure that Z (T^) in (23) has. The still obey the same 
balance condition = 0 as earlier. The total momentum is just as above Altotai = ~ mod Np. 

The generalized weight is now 


(r«) 


= exp < ZTTT 



nrTi 


ij z. 


{Kij) 


}> 


(51) 


which has the same type of structure as (23). Again the Tij are anti-symmetric just as for the Laughlin 
state. Similarly the factor is responsible of for the anti-symmetry under the exchange of two from 

the same group. 

Note that, as the different groups of electrons now are distinguishable, Z (T^) is not proportional to the 
Fock coefficients, which was the case for the Laughlin state. In fact the states in the hierarchy needs to go 
through (48) to become physical electronic wave functions for a single component state. As a result, the 


coefficients given by Z (T^) would not even be the Fock coefficients, but rather the components that enter 
into them. 


6.2. TT-limit scaling for HH-states 

In Section the asymptotic scaling of the Fock coefficients where extracted for the Laughlin state. A similar 
analysis should be possible also for the general hierarchy states using (51). Just as with the Laughlin state 


the dominant state should be the TT-state [TT83llBK08| . This is the state where the electrons are maximally 
separated along the long dimension of the torus. The analysis is made more complicated by the fact that 
Z is not invariant under Tj —> -|- Np which gives more than one weight that corresponds to the root 

partition. 

We can however, without performing that full analysis, shed some light on an observation made in an 
earlier paper. In Ref. |FHS14| it was noted that in comparison to the Coulomb ground state v = |, 
some of the translation operators that appear in (48) gives better overlaps than others. For instance, in the 
TT-limit, the translation operator Di q gives the best overlap with coulomb, whereas Dq i yields an overlap 
that is almost zero. See Fig. 3 of that paper. The understanding we had was that Di q constituted a smaller 
translation and therefore resembled a derivative more than Dq i did. From a technical point of view we can 
now also say that since the TT-limit is dominated by the single slater determinant of the root partition the 
effect of Di q and Hop are drastically different. All Fock-states are eigenstates of the operator Hi^oi and 
so it will only give rise to a phase factor. The Hop operator on the other hand will shift all the electrons 


in one group relative to all the others, which changes the configuration. 


As an example, for v = the 
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TT-configuration is... 010100101001010 .... When acted on by Hi o it will just pick up a phase and becomes 

7277 ^ . ^7 

e 010100101001010 .... Under the action of Dq i however, it becomes ... 001100011000110 ... 

which is a configuration with almost no overlap with the Coulomb ground state. 

7. Summary 

In this paper we have rewritten the Laughlin state in the momentum Fock-basis. The techniques developed 
here can also be applied to the full chiral Haldane-Halperin hierarchy wave functions on the torus, although 
the results are not as easily interpreted as for the Laughlin wave function. We have further re-expressed the 
Fock-coefficients of the Laughlin state in a recursive form for faster numerical evaluation. 

We used the analytic expressions for the Fock coefficients to study the Laughlin state deep in the TT- 
limit. We did this by analyzing the asymptotic r-dependence of the Fock coefficients. Using this asymptotic 
expansion we have analytically computed viscosity in the TT-limit as well as the leading r-dependence on the 
correction to the normalization. We found a correction to the proposed CFT-normalization that is accurate 
up to a factor f^‘ (L). 

Our analytical approximation for the viscosity show qualitative but not quantitative agreement with the 
numerically obtained viscosity. We conclude that the parameter that dictates the viscosity in the TT-limit 
is the short length scale L and not r. Thus, any state sufficiently deep into the TT-limit will always have 
the “trivial” viscosity of s = i. Or findings suggests that the plasma is not screening in the TT-limt. This is 
supported by the need to correct the normalization by a T 2 dependent factor when L < 4.4£b- This results 
thus warns that Hall viscosity can not be used as a probe to detect topological information on this asymmetric 
geometry. We still believe, however, that the plasma analogy holds in the thermodynamic limit, both for the 
Laughlin state and for the rest of the hierarchy states, and we have given some numerical evidence that this 
is the case in this paper. We further believe that you can study the viscosity as a function of L to determine 
whether a topological phase has stabilized or not. By thermodynamic limit, we here mean the limit where 
both the torus axes are much larger than the magnetic length. 

One of the advantages of studying the TT-limit is that it can be exactly solved while still being 
adiabatically connected to the bulk state at r = ijBKOH IBK08] . The gap to excitations remains open 
when tuning r between r = ^ and r = zoo and the charges of the quasi-particles are the same. What we 
have seen however in this paper is that some other properties does change, such as viscosity, but there are 
other quantities that also get altered. One such is the exclusion statistic |Hal91| for the quasi-particles which 
changes |KKll] at about the same L as the viscosity. This serves as a reminder that the TT-physics is in 
many cases the same as in the bulk, but that there are notable exceptions that one needs to be weary of. 
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Appendix A. A recursive formula for 

In this appendix we list some of the properties of the t?-function structure factor . This object appears 
in the exponentiation of the generalized t^-function as 



{z\rf 

_ 'y ^ ^i7rT^{k+aN)^ ^i27T{k-\-aN){z-\-b) 



tgz 

We define as 



II 

OO 

E 

f J. 1 

+2 


= —oo 


where ^^ U = 0 and ^ and ti = ii + ^. This enables us to rewrite it as 

OO 


(A.l) 


(A.2) 



Success and failure of the plasma analogy 


for Laughlin states on a torus 


28 


where letting U —>■ —ti we easily see that = Z^\ As one of the ki is linearly 

dependent on the others, we may without loss of generality we choose it to be tAr = fc — such that 

^ — In = ^ — k + ~ w ~ rewrite Z^^^ as depending on iV — 1 

linearly independent variables {ti}^ ^ 

OO 

where we introduced the dummy variable Tat-i = note that Z^^'^ = shift k —>■ k+N 

is reverted by the re-summation ti^ff —> t^^j^ + 1, i.e. Tat-i —> Tat-i + iV — 1. We now note that we may 
extract a reduced form of ( |W2| by identifying k with Tat-i and N with A^ — 1. We then get 

k / ^ T’jv-i 

tiV — 1 — — C)0 

Since Tat-i —>■ Tat_i + iV — 1 and k ^ k + N cancel each other we can rewrite the recursion using i?-functions 
as 


7{N) 


= E 

—1 


W-1)„ 


i /c * 

(o|7V(7V- 1)t) 


with the root function zj:^'^ = 1. 


(A.3) 
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Appendix B. Expanding Laughlin in a Fock Basis 

In this section we give the details on how to rewrite the real space Laughlin wave function (main-text:7) in 
the single particle orbital basis of (main-text:5), i.e. the Fock basis and arrive at equation (main-text:23). 

The key to constructing the Fock expansion is to extract the Fourier components from (main- 

text:?) as they will directly map onto the Fourier components of (main-text:!). In order to do so, these 
components need to be extracted from the Jastrow factors as well as the center of mass factor. At this stage 
of the construction one may identify all the factors containing x with a corresponding ()/j (z). Reconstructing 
the full single particle orbitals rjk {z) is then done by massaging the sums over k. 


Appendix B.l. Fourier expanding the Laughlin state 


In this part we give the steps to reach (main-text:23). When Fourier expanding the Jastrow factors we 
will use an expansion where the Fourier factor is explicit. The generalized r?-function to the power M can 
be expanded as in (main-text:25) where is a factor that encodes all the information about the 

power. We have the relations = Z^'^ and by mirror symmetry, {z\t) = 1)3 {—z\ t ) f function also 

y{M) _ y{M) 

Zj - Zj _ rp , 

Applying the expansion (main-text:25) to the full Jastrow factor in (main-text:?) gives the form 


n^i 

i<j 


{Tij GZ} 


|(T. + f)%*2^E.<,(T.+ §)(¥ + 


i)nz 

i<j 


(<?) 
T- • 5 


(B.l) 


where with each pair z^ has been associated with a summation index Tij. As Zij is anti-symmetric in its 
indices Zij = —Zji we will require that Ty be antisymmetric too. We choose a skewed anti-symmetry condition 
Tij + Tji = —q to make the anti-symmetry for Tij trivial. To facilitate further calculations we introduce the 
non-integer variable T^ = Tij -|- which has symmetry properties Tij = —Tji. From this follows that Ta = 0 
for all We also introduce the shifted Z-weight Z^\ = Z^?\ ^ = z!^^_ which has the same symmetries as 

Z^\ and especially obeys Z^\ = zjd^. We also introduce the variable 

= (B.2) 

i=i 

to simplify the z-dependent factor ^i^j TijZij = "Yh^j Zifij = Yf" Here is thus related to the total 
momentum of particle i. Note that because of the anti-symmetry of Ty there is a balance condition on 
giving 

We We 

^T,= ^r,,=0. (B.3) 

i—1 i^j — \ 

Putting all of these things together, the Jastrow factor can be written as 


i<j 


E ' 


Ei<,- Tij Y.i Ti ^ 


nz 


( 9 ) 

fa' 


i<j 


(B.4) 


t Here {z\r) = '& 


{z\ r) is the third Jacobi theta function. 


X Formally Ta is never defined, as za does not appear in the Jastrow factor i B.l I. 
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Note that we have managed to extract the phase factor ^ for each coordinate separately. The price 

we have paid is the introduction of the Tij which label the interdependence of the different momentum 


components. The balance condition (B.31 ensures that the wave function will be in a well defined momentum 
sector. 

We may in a similar fashion as for the Jastrow factor also express the CoM function in the relevant 
Fourier components. To be able to encompass any type of boundary conditions we work with the more 


general form iFh,t ({z} ,t) = d 


h 

Q 

t 


{qZ\ qr). At the end of the calculation we can set h = t = s and recover 


the result in (main-text:22). 

Using this parametrization the CoM piece is rewritten as 






(B.5) 


Putting together the expansions (|B.4[), 


i<j 


{f.,GZ+f} 

X ^ g*7rrg(m+|)^g227i-(Ti+m9+/i) 52; fr ^i2Tr(^mt+^f) ^ 


and the Gaussian factor gives us the wave function 

rj^ 2 


17TT - 

e ^^<3 q e 


(B.6) 


The next step is to extract the single particle wave functions (main-text:!) from (B.6) by forming basis state 
factors of the form (main-text:3). Me may now collect the factors containing Zi and yf into [zi) using 
(main-text:3). This gives rise to the substitution 


^ = e 






n 


i=l 


where 


ki = Ti -\- mq -|- h. 


(B.7) 


We can now see that the role that h plays is to ensure that ki is an integer. The parameter t plays a 
similar role but for the momentum in the y-direction. For periodic boundary conditions then ki € Z, which 
means that -|- h must be an integer. Since the offset on is f “ ~ 1) § term T), = 0 does 

not contribute) the same offset has to apply for h. Expanding the counter weight term gives 


—l-KT - (mq+h)^ 

e ’ = g "0 » g g\ 'i r , 


where we used (B.3l. Putting these terms back gives 


{UjGZ+f} 

X v) n 4^) n a. (^^) • (B.8) 

i<j ^ i=l 


Before we proceed let us stop and inspect the current state of 'iph,t- The expression (B .81 already has 
some of the desired structure. For instance we have managed to isolate the factors fki (zi) which are the 
building blocks of the basis states in (main-text:4). In order to get to (main-text:22) though, we need to 
perform parts of the sums over Tij. This is what we will do in the next section. 


Appendix B.2. Fock coefficient for the Laughlin state 

Here we give the extra steps needed to arrive at the Laughlin state (main-text:23). The remaining task is to 
make ki independent of m such that ki = Ti + h. 

Note that for periodic boundary conditions we mush choose t and h as the half-integer t = q {Ng — 1) | 
and h = q {Ng — 1) ^ + Z. This will work for both fermions and bosons and in the latter case h and t would 
always be integers. 
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We start by grouping (B.8l as 




'^h,t = 






where 
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T2 
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(B.9) 


(B.IO) 


i<j 


Here we have dropped an overall factor of 9 . This is the same Z (T) that appears in ( main -text:23). 


and the Tij are still subject to the constraint (B.2|. By extracting Z (T) we have formulated (B.9l in a very 


suggestive form reminiscent of (main-text:22), only that the summation is a little off. For instance have the 
extra factor of a sum over m and the sum over is infinite. We remove by shifting all the 

Ti by Ti —)■ Ti — qm except for which by the constraint (B.3l goes to —>■ + (W — 1) qm. This 

is performed in practice by shifting the sums over the as T^- —>■ Tij + q'm{5i^i<i^ — Sj^N^) ■ Under this 
transformation the momenta change as 

ki —y 'Ti h 

and Z (T) —>■ Z (T) This extra phase extracted will exactly cancel the factor present in 

(|B.9p. We also split as Tj —>■ Ti+riN^ where now G and € Z, such that also obeys the balance 
condition J^i = 0- U can be shown that Z (T^ -f riN^) = Z (T^) is independent of r^. Finally shifting 
m ^ m — enables us to complete the sum (main-text:4) to obtain the expression (main-text:22), where 
ki = Ti + hmodN^. It can be shown, as a sanity check, that Z (T) is fully anti-symmetric in th interchange 
of its variables Ti. 


Appendix C. A Recursion Formula for Z 


In section Appendix B we showed that the Laughlin state (main-text:7)could be expanded in a Fock basis 
with coefficients Z (T) given by (main-text:23). Unfortunately the given expression is not particularly helpful 
when it comes to numerical evaluation. The culprit is the simultaneous infinite sums over all of the Tj^cs. To 
alleviate this problem a bit we will in this section rewrite Z (T) such that the sums can be performed in a 
recursive manner. This will substantially reduce the scaling of the computation. 


This section is divided into three parts. The first part. Appendix C.l deals with rewriting Z (T) for A, 


particles such that it depends recursively on the components for — 1 particles. The remaining two parts, 
further manipulates the expression for Z (T) for enhanced numerical efficiency. 

A word of caution should be given. Although the final expression (C.lOl can be put on a computer. 


the efficiency in evaluating it is still inferior to that of simply diagonalizing the Haldane pseudo-potential 
Hamiltonian. Using exact diagonalization the coefficients for around 12 particles can be extracted numerically. 


but in our current implementation of (C.IO), only 6 or so particles can be achieved. The reader looking at 


equations (C.IO) and (C.lll may be worried that they are wrong or will render errors when implemented 


numerically. We have however validated our code against the exactly diagonalized Haldane pseudo-potential 
Hamiltonian and have perfect agreement for all tested cases. 

With the given word of caution in mind, the reader should also be aware that the pseudo-potential trick 
only exists for the Laughlin state. For states higher in the hierarchy there is no local Hamiltonian for which 
these state are the exact zero energy eigenstates. In these cases the results in this article is the only way we 
are aware of to analytically extract Fock coefficients. 


Appendix C.l. Recursion in 

In this part, we write 2 (T) in a recursive formulation such that the coefficients for N particles depends on 
the coefficients for A — 1 particles. For analytical purposes the Gaussian exponential piece of Z (T) can be 
written in different ways as 


i i<j<.k 


i<j 
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i<j 


T, - T, 
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Especially the two last expression shows explicitly that Z (T) is a converging sum, but they do not 
particularly help when numerically evaluating Z (T). In order to find a numerically tractable formulation of 
Z (T) we will develop a recursive construction in N. For this purpose we introduce some more notation. We 


write zi^'^ as the Fock coefficients for N electrons. We also write 


- T 


to keep track of which 


electrons are being included in the momentum at level N. This induces the relation between different 
system sizes as 


f.AT (C.l) 

By extracting all the pieces from Zj^'^ that belong in Zj^ it is possible to recursively formulate the former 
in terms of the latter. This step is straight forward as any factors containing T^- where both i,j < N can 
be moved to Zj^ immediately. The remaining pieces are TiN which can be rewritten according to (C.l I 
such that the recursive form of zi^'^ is 


r(N) 


= S{N) 
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^z7rrq'A(A^) gZTT ^ 1) 


N-1 


.{N-l) 


GZ+g- 


i} 


n ?'■" 




(C.2) 


In the above equation we have introduced the difference between the Gaussian exponentials as 


11 N-i / 

A(N) = ^A{N)-^AiN-l) = iN-l)Nj2 

9 1 i=i \ 

The prefactor S {N) = 
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(Af-l) 


t: 


(TV) 
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(TV)2 
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qiN-1) q{N)J q^N ' 
T is the quotient between the two factors ^ 


(C.3) 


and 


(AT-l)^ 

AT-l 

^ 4 Already this re-formulation has reduced the computational complexity of Z^^'> from order 

O since each pair of indexes had a sum Ty, to a still hard but more humble O (e^)- This is of 

course provided that the results further down the recursion can be stored and reused. To further facilitate 
the evaluation we seek a formulation that as much as possible disentangles the different sums from 

each other. This will the purpose of the next part. 

Appendix C. 2. Orthogonalizing the sum in A (iV) 


With the form of A {N) in (C.31 we can compute recursively. There are two computational problems 

that need to be remedied though. The first is the form of A (N), as the different sums over T^-'^ are not 
separable, i.e. they contain cross terms T^Tj. These cross terms reduce the numerical efficiently as all the 
sums need to be evaluated simultaneously. The second problem is that the sums over are infinite 

whereas Z^^ are defined modulo q {N — 1). We thus seek to rewrite A (N) such that this modularity is 
explicit. In what follows we will remedy both of these problems in turn. 






To reduce the notational complexity we define = Ti and = Ri such that (C.31 can be written 


as 


^ (^) ^ (r - R 

N{N-1) A N_i 


■{N) = 


(C.4) 


(TV) 


are subject to neutrality conditions Ri = 0 and ^ = 0 and 


Note that Ri and Ti - just as 
herein lies the difficulty. The neutrality condition means that (C.41 is not diagonal in neither Ri or T^. The 


first order of business is to rewrite (C.41 on an explicitly positive definite form. The details of the rewriting 
are in [Appendix E| and the result is 


M 

TT {M -I- 2) = w„ 

n—1 
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C’l) - V E 
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We chose to use M here instead of — 2 to empathize that there are M independent variables Tj. The 
vectors are 


(**) / ■ I 1 

Uj = \ —n J = n + l 


1 


j <n 


(C.6) 

0 j > n + 1 

and are eigenvectors to the M x M matrix A4ij = 1 + 5ij. The matrix has eigenvalues A„ = 1 + M6n^M- 
The squared norm of the vectors is = n + n^ (1 — 5 m, n) which enters (C.5l as Wn = ' 


As (C.5l stands, the sums over Tj are still entangled. Note that there exists no integer shifts of the Ti 
that will diagonalize the sum. The simple argument for this is because det = Ml, which means that 

the inverse of is not integer valued. We note however, that the sum be manipulated 

to the form (1 + n) T„ — nTn+i (n < M) or Tm {n = M) by the simultaneous shifts Tj —)■ Tj — Tj_i for 
1 < J < M. We will use this fact in the next part to reduce the complexity in evaluating . 


Appendix C.3. From tt {M + 2) to d-functions 

In order to arrive at an efficient computation we must make contact with the periodicity of {N — 1) q in 
Zj^ This periodicity coincides with —>■ +Z. For this purpose we therefore extract exactly this piece 

from Ti by splitting of the fractional pieces as Ti —)■ Ti 


Here T, e Z and e Z 


g(M+l) 


qM 

2 ■ 


We proceed by simultaneously shifting Tj —>■ Tj — Tj-i for 1 < j < M. The positive definite squares are now 
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■(M + 2) = ^ 


Wn 
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(n + 1) r„ - nTn+i - + WM (tm - 

'^n J V 


where Dn is the combination of Ri and giving 


Dn=-Y. 


(n) 


^ ^ q{M+l) 


M+1 
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(n)_ 

q{M + 2)' 


(C.7) 


(C.8) 


We further manipulate (C.71 by separating T„ as Tn —>■ nln + qn, where In & and qn G Z„. For the term 
with n = 1 this trivially becomes Ti —> /i. On top of this we remove the cross terms between In and In+i 
by shifting In ^ In + ^n+i starting from n = 1 and going upward to n = M — 1. More details are given in 


Appendix F The expression for tt (M + 2) now becomes 


M j n \ ■ 

(M + 2) = \2{n+l)n(ln + --j^^- , 

V « (n+l) n(n + l); 


(C.9) 


where we have introduced the dummy index qM+i = 0 to write the whole expression as one sum. The 
sums over In can now be performed to produce il-functions, such that the full recursive product is 




= S{N) 




X e 
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The introduced weight is 


N-2 N-2 




^n + l 

n+l 


Dr, 

i(n+l) 


(0| (iV - 1) iV (n + 1) nqr) .(C.ll) 


The il-factors in (C.ll I converge fast since exp (ztt {N — 1) N (n 1) nqr) will usually be a small number. 
Also the sum over J2{q~^z } looks naively as is contains a number of terms which will scale as {N — 2)!. 

However, because of the structure ^ It can be decomposed into n{n-\-1) = O {N^) evaluations 

of -d-functions. 
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Appendix D. Fock Expansion of the Hierarchy States 

In this part we give the steps to reach (main-text:72). The construction will parallel the construction for 
the Laughlin states albeit with more structure. We will here only point out the main differences and display 
important steps in the expansion. 


Appendix D.l. Fourier expanding the hierarehy 

Just as for the Laughlin state we begin with the Jastrow factors. Again we use (main-text:25)to expand the 
i?i-functions. The expansion reads 






i<j 






g - ICij E: 
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T’ ’ 


i<j 


(D.l) 


where just as before T^j = —Tji. Note that since the pair Zij now comes with exponent there are 
different types of Zt present. However, they still obey the same equations Zt 

and Note that has its maximum TT-value at Similarly as before, we have 

ii ii ii — 2-3- 


_ rzK^ij) _ 


also introduced the variable accordance with (B.2). The balance condition on is the 

same as earlier ~ The balance conditions (B.3l ensures that the wave function will be in a well 

defined momentum sector. 

We may in a similar fashion as for the Jastrow factor also express the CoM function in the relevant 
Fourier components. This part is more complicated compared to Laughlin case there is more structure in 
the CoM term (main-text:67). 

For this purpose we parametrize the charge lattice vector q as q = maq_a where {qa} span the charge 
lattice F. Since we are only considering the abelian chiral hierarchy it suffices to use the one-dimensional 
parameterization h = hho and t = tho. At the end of the calculation we can set h = t = s and recover the 
result in (main-text:22) and (main-text:72). Here ho = has the property ho • q^ = 1 and hg = ^ = |. 
This parametrization is generic enough to capture all the states in the abelian chiral hierarchy. Note that 
gho G F which means that all q degenerate states at v = ^ are obtained by letting s take values from 1 to 


q. For states with higher degeneracy than q, such as the Halperin 331-state with K = 


3 1 
1 3 


a more 


elaborate parametrization will be necessary. For more detailed properties of ho and F see e.g. Appendix B 
of Ref. |FHS14| . 

Using this parametrization the center off mass piece (main-text:67) can after expansion be written as 


{me} 


Again, when we combine the expansions (D.l I, (D.21 and the Gaussian we may isolate the single particle 
wave function building blocks (main-text:3). We substitute for fkt {zi), but this time there is more 

structure in the momentum relation 


ki =Ti + maKai + h. (D.3) 

Oi 

As earlier we see that both h ensures that ki is an integer. For periodic boundary conditions then ki S Z, 
which means that Ti + h must be an integer. Since the offset on is —~ term Ta = 0 

does not contribute) the same offset has to apply for h. We are here assuming that all Kaa have the same 
parity, a reasonable assumption equivalent to demanding that either all particles are bosons or all particles 
are fermions. 
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Substituting for fki (zi) and expanding the counter weight term e fcf using 

the expansion reads 

X [| Cfc, (z*). (D.4) 

i<j i=l 

In this expression we now have extracted the Cfei for the individual particles. Not also that h no longer 
is a parameter in the expansion. This is to be expected as the only difference between the q different states 
is the labeling of the Fock coefficients. 

Appendix D.2. Resummation of Tij 

We now continue and perform some manipulations on the sums of Tij. The aim of these manipulations is to 
collect the rua in ki on specific particles. We note that ki = such that the could 

be used to complete the construction of the basis functions rjk {z). It can be shown that the shifts 


B.31 and Q-qj = Nj, 


^ (5, - Sj) (D.5) 

rric^ ^ Si^ , 

iaGla 


will send ki ^ ki + Nj,Si while leaving the rest of (D.41 invariant. This shows that it is in principle possible 
to construct the rjk (z) for all the particles by shifting the sums in appropriate ways. 

We will however here choose a slightly different shift on the summation. For this purpose we shift the 
Tij as 


^ij ^ '^ij T ^ ^ ' 

a 


This transformation only affects the indices for the “last” particle in each group, ^411 and has the effect 
that ki ^ ki = Ti + Nj, rriaSi^Na + Iii effect it puts all the ma on only one particle per group. This 
transformation does nothing to as it is periodic under Tij —> Tij + Kij. 

An extra phase of exp (itt rua (Nj, — Kaa)) is picked up which precisely cancels the existing phase 


exp (zTT mat). There are extra factors picked up from 
come together in such a ways are to precisely cancel all 
things together we have the expression 


'Hici 


and 


which 


dependence in the Fock weight. Taking these 


i’h.t = X! ® 


-tj g 


Ne 




r(K^, 

1IJ 

'^<j {rria} i—1 




(D.6) 


where we have also dropped the constant phase e . Note that as ttIq, is not part of the weight, the sums 

over TOq can be performed to produce n of the W rjk {z) needed to construct the Fock coefficients. In coming 
this far we can now identify the factor making up Z in (main-text:72), and the rest is (main-text:70)given in 
the main text. 

t Note here that we abuse notation and let Na be the index of the last particle in group a. instead of the size of group ol. 
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Appendix E. Prom (C.4) to (C.5) 


In this appendix we give a detailed derivation leading to equation (C.51 in the main text. We start from 


equation (C.4) in the main text: 

N-l 


R 




Z=1 


(E.l) 


Note that Ri and R also are subject to neutrality conditions Rj = 0 and ^ R, which means that 


^N-l , 


(E.l) is not diagonal in neither Ri or R. We begin by rewriting ( |E.1[ ) on an explicitly positive definite form 
with only N — 2 terms. Since the number of linearly independent variables R variables are iV — 1 we also 
define M = N — 2 to keep track of this number. 

We proceed in two steps. First we set set Ri = 0 and diagonalize only T^. After that we insert Ri again 
and deduce the full expression. 

Appendix E.l. Simplification Ri = 0 
We start by defining the simplified 

M+l M / M \ ^ 

7To{M + 2) = 7T{M + 2)\^^^,= Y,Ti=Y.Ti+iY.^A ’ (E.2) 

\i^l / 

where only the Ti are present. In matrix form this would correspond to 

M 

^0 (M + 2) = ^ (E.3) 

* j'=i 

with the matrix Mij = 1 + Sij. The eigenvalues of that matrix are A„ = 1 + MSm.u- The eigenvectors can be 
constructed by noting that = 1 is the eigenvector to Am with eigenvalue M + l. The rest of the vectors 

form an degenerate space where all vectors have eigenvalue 1. The full space of eigenvectors can therefore be 

parametrized as 

1 j <n 

= —n j = n + l . (E.4) 

[ 0 j > n + 1 

The squared norm of these states are 

s„ = n + (1 - Sm.u) ■ 

In terms of these ttq (M + 2) can be rewritten as 




M 


M 


TTo (M + 2) = ^ u;„ I ^ 

n=l \i=l 

_ Xn _ + 


(E.5) 


where Wn = = -— n(n+i) ■ remind ourselves of the special case M = 1 in which ttq (M + 2) = 2T{. 

Appendix E. 2. Restoring the factors of Ri 

We now restore the parameters Ri. We can not however put them back by simply shifting R ^ R — Ri, as 
the term would be handled in the wrong way. Instead, we deduce the values of these terms by 

studying the cross terms TiRj. Let’s take a look again at tt (M + 2) in (E.l) but this time writing out all the 
linearly independent terms of both R and Ri: 

M / M \ 


■{M + 2) = Y,{T,-R,f+[Y,T^ + RM+l] - 




R: 


M + l 


(E.6) 


To match this we make the ansatz that all of the Ri factors can be accommodated in the positive definite 
squares as 


M 


M 


■ (M + 2) = ^ I ^ - Ar, 


(E.7) 
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We here assume that 

M+l 


An =Y. ’ 


i=i 


is some linear combination of Rj. By compa ring the two expressions for (E.6l and (E.7l we deduce the 
coefficients By expanding the squares in (E.6) we have 


M 


M +1 f ^ ^ 

7^{M + 2) = 7^oiM + 2) + 2J2T^ {Rm+i - R^) + Y. ^ 


M. + 1 
i~l 


2=1 


2 = 1 


Doing the same expansion for (E.7) gives 


MM M 

,(«) 


' (M + 2) — TTq {M + 2) — 2 'y ^ WnAn ^ ^ Tj + ^ ^ Wn-A^. 

n=l j=l 


(E.8) 


(E.9) 


n=l 


As the two expressions for tt (M + 2) should be equal irrespective of the values of R we get equations for A„ 
that are 

M 


Y. Anvf'^Wn 

= Ri — Rm+1 

(E.IO) 

M 



y , '^n^n 

n—1 

n2 / 

* M + l 

2 = 1 

(E.ll) 


The upper line comes from the linear terms in R and the lower equation is the constant piece. Equation 


(E.lOl is really a system system of equations M unknowns and M equations. The solutions are 


.. M+l 

An = 

An . _ 


Then, putting An back into the ansatz we arrive at we expression 

2 

M I M ^ M+l 


■ {^M + 2) = Y^n\Y^iA-Y:ll^i^^\ , 

n=l Vf=l 


i=l 


with is expression (C.5l in the main text. 


Appendix F. Prom (C.7) to (C.9) 


In this section we give the detail of how to reach equation (C.91 from (C.7l. We may rewrite it slightly such 
that it is 


1 

^ (M + 2) = Y . , ((n + 1) r„ - nT„+i(5„^M - . (F.l) 

In the expression above we have abused notation and wrote Sn^M = 1 ~ ^n,M- Note that each squared term 
only contains two R terms, but that they come with different prefactors. As the context of this equation in 
insensitive to integer changes in R, we may shift the summation at will without needing to consider other 
external factors. We proceed by splitting sums over r„ in two pieces atr„ = nR + Qn where G Z„ and 
In € Z. For n = 1 this trivially means R = li. Doing this results in 

M 

7t{M + 2) = Y ! 

^nin+l) 

n—1 ^ ' 


((yr “h 1) U (^n ^n+l^n^M) T ((^ T 1) Qn ^Q^n+l^n^M) ^n) 
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The facts that are summed to infinity, i.e. In now come with the same prefactor n (n + 1). We can thus 
eliminate one of the terms by shifting In ^ In + In+i starting from n = 1 and going upward to n = M — 1 . 
This will yield the expression 


M 

TT {M + 2) = (n + 1) n 

n—1 



gn+l 
n + 1 




Dr, 


. (n + 1) 


2 


By treating qm+i as a dummy index qu+i = 0 we can remove the Sn^M and obtain the homogeneous looking 
expression 


M , 

TT (M) = {n + 1) n (In + — — 
V n 

r> = 1 ' 


gn+l 

(n+ 1) 


which is (C.9) in the main text. 


Dn y 
n (n + 1)/ 


(F.2) 







